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We study the dynamics of a pair of molecular ensembles trapped inside a superconducting res- 
onator through which they are strongly coupled via a microwave field mode. We find that entan- 
glement can be generated via "vacuum fluctuations" even when the molecules and cavity field are 
initially prepared in their ground state. This entanglement is created in a relatively short time and 
without the need for further manipulation of the system. It does, therefore, provide a convenient 
scheme to entangle two mesoscopic systems, and may well be useful quantum information processing. 

PACS numbers: 



Vacuum fluctuations can have important physical con- 
sequences, for example, in the Casimir effect [l| and 
Hawking radiation [2|. van der Waals interactions, i.e., 
attractive long-ranged interaction between neutral atoms 
or molecules, are also a kind of Casimir effect. It is an 
interesting question as to how vacuum fluctuations can 
be used to influence the properties of quantum entangle- 
ment between two systems. Quantum entanglement is 
a fundamental concept in quantum mechanics Q and is 
also thephysical resource in quantum information pro- 
cessing 4]. In fact, it has recently been shown possible 
to generate entanglement between a pair of particles via 
the vacuum modes of the radiation field [!, E| 0] • 

In this paper, we study how vacuum fluctuations in- 
duce quantum entanglement between two mesoscopic sys- 
tems, i.e., polar molecular ensembles Q are placed in- 
side a cavity, and strongly coupled by a single microwave 
mode. Recently, Rabl et al. [§] have proposed the re- 
alization of a quantum memory using such ensembles of 
polar molecules inside a superconducting resonator [lfl ]. 
The energy difference between two internal states of a 
polar molecule is the order of GHz and polar molecules 
have significant electric dipole moments. A strong cou- 
pling to a microwave field via a transmission line can thus 
be achieved. In addition to the strength of the coupling, 
low-lying collective excitations can be coupled to the field 
and exploit the enhanced coupling to them, which scales 
as \/N, where N is the number of molecules in the en- 
semble. 

The dynamics of vacuum fluctuations [1 ll ] is hard to 
observe in ordinary systems. To show why this is so, 
we start the simple case of a two-level atom interact- 
ing with a quantized field. Conventionally, we use the 



Jaynes-Cummings model [12| in the interaction picture, 



Hamiltonian H = g'[tf a_e i{ - w '-"o)* + cr+fee" 4 ^'""^*], to 
describe a two-level system a± coupled to a quantized 
field 6, for lu' Qi uj 1 and g' are an energy difference between 
two-level atom, the frequency of the field and the Rabi 
frequency respectively. The rotating wave approximation 
(RWA) can usually be used because the two countering- 
rotating terms, 6cr_ and <J+b\ can be neglected; they 



carry a fast oscillation with the high frequency lu' + lu' q . 
The RWA is, therefore, an excellent approximation for 
the optical frequency regime in the weak Rabi coupling 
limit. Clearly, this Hamiltonian will produce no evolution 
in the atoms and the photon field if they both start in the 
irrespective ground states. However, this approximation 
breaks down if the Rabi frequency g' is comparable to the 
frequencies uj' and ui' . In fact, the RWA is completely 
inadequate to describe the physical situation of a large 
number of molecules interacting with a microwave field 
in the strong coupling regime. It is thus necessary to go 
beyond RWA and, in essence, study the role of vacuum 
modes on the dynamics of the coupled atom-field system. 
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FIG. 1: Two molecular ensembles, separate with a distance d, 
are trapped inside a superconducting resonator and interact 
with a single-mode microwave field along the stripline. 

We consider the case where molecules and the pho- 
ton field arc initially prepared in their ground states and 
show how the countering-rotating terms in the Hamilto- 
nian do indeed induce quantum correlations between the 
molecules. The two ensembles of molecules exchange en- 
ergy with the vacuum field due to the counter rotating 
terms. In physical terms we would say that this proceeds 
via virtual excitations of the cavity mode. We should 
bear in mind that these vacuum mode processes can be 
enhanced if collective excitations of the ensembles are 
used. As the dynamics takes place, the two molecular 
ensembles becomes entangled as the molecules are effec- 
tively coupled through the exchange of energy via the 
vacuum mode of the cavity. We will show that this en- 
tanglement can be generated in a comparatively short 
time. This result provides a novel route to entangle two 
mesoscopic systems. Relaxation and decoherence effects 
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should also be modest as the system is prepared in a vac- 
uum state; a potentially crucial advantage for practical 
applications. 

In this paper, we suppose the molecular ensembles 
are placed at a distance d apart in an effectively one- 
dimensional resonator as shown in Fig. Q] The two states, 
produced by the splitting of a rotational ground level by 
hyperfine interactions are chosen as the internal molec- 
ular states ||. They have an energy difference u, and 
are coupled to a single-mode microwave field with a fre- 
quency ujq. Each of the molecule can be described by a 
spin-half particle <Ji, and hence a collection of such spin- 
half particles can be described by a collective angular 
momentum operator J = ®i-> w hen N is the num- 

ber of molecules in one of the ensembles. The wavelength 
of microwave radiation will be much longer than the size 
of molecular ensemble. Hence, we can then assume the 
microwave field coupling to all molecules for the system 
of coupled molecules and radiation field in the form. We 
can now write down the Hamiltonian H (h = 1) [l3j |: 

2 

H = J2LO a f a + LuJ Zi + c H {a + a))(J +z + J.J, (1) 
i=i 

Here, and a is the creation and annihilator operators of 
the cavity mode, J Zi and J Xi are the angular momentum 
operators to describe the collective inversion and tran- 
sition for the i-th ensemble respectively, and i = 1,2. 
The molecule-photon interaction strength is denoted by 
gi for the z-th ensemble and they differ with a relative 
phase <p = 2nujd/c between the field and two ensembles, 
where c is the speed of the microwave photon. For sim- 
plicity, the magnitude for two Rabi coupling strengths 
are chosen to be the same, \gi\ = \g2\, and </>«0. We con- 
sider the case where the molecules and photon field are 
in resonance, i.e. uj = luq, which is the optimal condition 
to observe the effect of the small vacuum fluctuations. 
We note that the Hamiltonian H has the same form as 
the Hamiltonian of the Dicke model without the rotat- 
ing wave approximation 14 1. The analysis we present 



here applies to thermal atomic ensembles as well as con- 
densates [15J. Condensates would have the advantage of 



longer coherence times but also introduce nonlinear dy- 
namical problems [l6| ■ 

To study the quantum dynamics of this system, we 
need to find the eigen-states of the whole system. The 
Hilbert space of this system is extremely large as the 
number of molecules are correspondingly large. We 
can, however, consider the time evolution case where in- 
volves only the low- lying excitations of the molecules. 
We can make this approximation based on the Holstein- 
Primakoff transformation (HPT) 17[ which enables us 
to express angular momentum operators in terms of har- 
monic oscillators. In this manner, we can take the leading 
approximation and map an angular momentum operator 
into a harmonic oscillators by taking the lowest order ver- 
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FIG. 2: Time evolution (the dimensionless time gt) of the 
logarithmic negativity lniV with the number of atoms N being 
10 4 . The different cases of uj = 300#, 500# and 2000g are 
shown in solid, dashed and dash-dotted lines respectively. 



sion of the HPT [3. We then get J Xl w ^ftJN~x a jV2, 
J Vi w -VNlPajV^, J Zi w (p 2 ai +u 2 xI z )/2uj-N 1 I2 for 
x ai and p ai are the position and momentum operators, 
and Ni is the number of molecules in the cavity. This ap- 
proximation is valid as long as (p 2 . + u 2 x^. (t)) / 2u;Ni 
[18|. It should be a very good approximation for the 
number of molecules are sufficiently large [9( . 

For convenience, we write the cavity field operators in 
the phase-space representation: x c — (a^ + a)/ \/2ujq and 
p c = i^/UJoi^ ~ a )l v2- We represent the system in terms 
of position and momentum operators. The Hamiltonian 
of system can then be rewritten in the form: 



H' = ±£te 



2 +pI+UoxI+u 2 xI+< 



i\J NiLJU Xa t X c ). 

(2) 

We now want to find the dynamics as the ensembles in- 
teract and become entangled. This problem is clearly 
related to what of finding the entanglement between two 
harmonic oscillators in an open-ended harmonic chain 

This harmonic system will be in a Gaussian state, 
which allows us to quantify the general bipartite entan- 
glement of the system. The density matrix of a Gaussian 
state can be completely determined by the second-order 
moments of the position and momentum operators of the 
system. We just need to study the reduced density ma- 
trix of the molecular part to find out the entanglement 
produced between the two ensembles. This reduced den- 
sity matrix can be obtained by the tracing out the cavity 
mode. The reduced density matrix p ai2 with matrix ele- 
ment (XiXj + XjXi) —2(Xi) (Xj) , where Xj — x ai orp ai - 
A quantitative measure of entanglement can be obtained 
by using the logarithmic negativity [2(| which gives us 
an upper bound for distillable entanglement. The loga- 
rithmic negativity in a Gaussian state can be found as 
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InJV = -5^1og 2 [min(l, 



7; I 



(3) 



where are the symplectic eigenvalues of the matrix 

We are now ready to investigate the entanglement dy- 
namics of this system. We consider the initial state as 
the state of molecules and cavity, i.e., the state of the 
decoupled harmonic oscillators. In Fig. [2j we plot the 
time evolution of the entanglement between the ensem- 
bles. The system begins with a separable states and then 
the entanglement grows rapidly. In fact, the quantum 
state of two ensembles oscillates between being separable 
and entangled. This is very similar to the entanglement 
propagation in a harmonic chain in which the two oscil- 
lators are periodically entangled and unentangled [HI . 

Moreover, the system achieves the first maximal en- 
tanglement within the time interval t* = 5~ 3 <? _1 . We 
can estimate this time t* with the realistic parameters. 
If we take g as 1 MHz nearly maximal entanglement 
can be yielded within 5 ns. This means that a significant 
entanglement can be obtained rather quickly. Moreover, 
no further adjustment of the experimental parameters or 
making conditional measurements are required (2l| . The 
time scale of this entanglement generation (~ 1 ns) is 
much shorter than the other decoherence sources such as 
inelastic collisions (~ 10 /xs) and photon losses (~ 1 fis) 
0. Entanglement can therefore be observed before de- 
coherence effect set in showing in a natural and efficient 
way to generate quantum entanglement for two meso- 
scopic systems. In addition, we can see that a larger 
ratio of g/u> can produce a larger degree of entanglement 
in Fig. [2J clearly indicating that counter-rotating terms 
cannot be neglected in this strong coupling limit. 

We should note that thermal noise is the main po- 
tentially problem in entanglement production. It is of 
course impossible to prepare the perfect vacuum state of 
a molecular ensemble in an experiment due to finite tem- 
perature effects. We now assume these ensembles can be 
described as a thermal state with mean excitation num- 
ber n — [exp (hiv/hsT) — l] -1 , for ks is the Boltzmann 
constant and T is the temperature. We can estimate n 
to be of order 0.1 to 1 when io ~1 GHz and T ~ 1 to 
10 mK 9(. From this estimation, we can see that ther- 
mal effects cannot be neglected, and it is important to 
study their influence on entanglement. Time evolution of 
entanglement under the thermal effects is shown in Fig. 
[3l The amount of entanglement produced is shown to be 
lesser in the cases of higher n. Moreover, the longer onset 
time of entanglement is required as shown in the higher 
temperature cases. But the entanglement can still be ob- 
served even if n is as high as 0.2. This result shows that 
a substantial amount of quantum entanglement can be 
effectively produced using thermal ensembles but colder 
molecules due result in a much better performance. 



In TV 




0.025 



FIG. 3: The logarithmic negativity In N is plotted against the 
dimensionless time gt with u = 300<? and TV = 10 4 . The solid, 
dash-dotted, dashed and dotted lines are represented n = 0, 
0.05, 0.1 and 0.2 respectively. 



Having discussed the production of entanglement, we 
now study how to observe the quantum correlations. In 
this Gaussian system, the density matrix can be con- 
structed if the uncertainties of these two ensembles can 
be obtained. This means that the entanglement of the 
two molecular ensembles can be determined just from 
the quantum uncertainties. In fact, non-resonant stim- 
ulated Raman scattering has been used to generate and 
verify the entanglement between two atomic ensembles 
[2 lL I22L [23j. In this scheme, the Stokes pulses are used 
to "write" the quantum information on the atomic en- 
sembles and the scattered Stokes photons carry the in- 
formation of excitations of each ensembles. Then, the 
two Stokes photon fields coming from each ensemble pass 
through a 50:50 beam splitter (BS) so that the two modes 
interfere and mix together. The conditional measurement 
of the resultant Stokes field can be pre formed and en- 



tangle the two atomic ensembles [2l|, [22[ . Similarly, the 



anti-Stokes pulses can be applied to read the excitations 
of the atoms and then the entanglement can be verified 
by measuring the correlations of photon fields. 

We now adopt this method to determine the entangle- 
ment of the two molecular ensembles as shown in Fig. 
[U We can apply two anti-Stokes pulses on these two en- 
sembles (being vertical to the view of Fig. [TJ) by passing 
through two small pinholes of the cavity. These would 
have to be larger when the wavelength of the optical ra- 
diation used for the probe and therefore should not affect 
the quality of the microwave cavity. To read out the exci- 
tations of the molecules, we apply an anti-Stokes pulse to 
each ensemble which optically pump the excited states to 
a higher rotational states with a large detuning between 
the ground rotational state. In the Hesienberg picture, 
the output beam, after passing through the ensemble , 
is given by af — Jrfia + y/1 — rfia{ , where af and 



the output field, the vacuum field operator and r\i is the 
effective transmission efficiency respectively. We can see 
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FIG. 4: In (a), Illustration for homedyne detection of measur- 
ing quantum correlations of a single molecular ensemble. Two 
anti-Stokes pump pulses are applied onto two molecular en- 
sembles. The output field from molecules is superimposed on 
the field from a local oscillator (LO) at a 50:50 BS and then 
they are detected by the detectors Di and D2 respectively. 
In (b), The output fields from two molecular ensembles are 
interfered with a 50:50 BS and then the fields are probed by 
the balanced homodyne detection. 



that the output field directly carries the information of 
collective excitations of the molecules. Hence, the quan- 
tum state of the ensembles can be determined through 
the measurement of photon fields. 

The measurement of the local and joint quantum cor- 
relations of molecules enables us to determine the log- 
arithmic negativity. We then require the detection of 
individual ensemble and the two ensembles respectively. 
In Fig IHa) , we give a sketch of a scheme to measure 
the local quantum correlations of an individual ensemble 
by balanced homodyne detection method via inputting 
a local oscillator mode of a coherent state with a large 
amplitude and phase </>j. [13]. The moments (xf), (pf) 
and (xiPi + piXi) can all be probed by appropriately ad- 
justing the phase angle <pi. Similarly, the joint quantum 
correlations can also probed by this method [12(. This 
can be done by interfering two output fields with a 50:50 
BS and then performing balanced homodyne detection as 
indicated in Fig. 0|b) • The quadrature of the two modes 
can be thus determined. 

In summary, we have found an efficient method to gen- 
erate entanglement between two separate ensembles of 
molecules and proposed a method to measure it. We 
have assessed the role of the finite temperature on to 
the entanglement produced. It is useful to the quan- 



tum information processing with molecular systems in 
a superconducting device. Our study has implication 
of quantum optics of mesoscopic system in the strong 
coupling limit. We envisage that evaporative cooling of 
the trapped molecules will realized [2J] so that the tem- 
perature can be lowered and the performance of quan- 
tum memory and entanglement generation be further im- 
proved. 
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